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Nomenclature 


© = Bearing radial clearance 
D = Fluid bearing damping 
D. = Rotor system structural and external damping 
F = Poree in + y direction 
Ka = Fluid film Direct Stiffness 
K, = Fluid film Quadrature Stiffness 
Ky = Rotor effective (modal) stiffness 
Mo = Rotor effective (modal) mass 
x = Distance in x direction 
y = Distance in y direction 
A = Fluid average circumferential velocity ratio 
+ = Shaft attitude angle 
$ = Angular velocity of rotor precession (Lomakin) 
& = Rotor rotative speed 
w = Angular velocity of rotor precession 
wa = Kotor system natural frequency 
W = Rotor weight 


here have been many workers in machinery vibration 

and rotor dynamics who have contributed to the 

advancement of the field. I can easily name several 
people in industry or in academia, in many countries of the 
world, that I look to as my senior in capability. 

Many of what seem to be “new” inventions are really devel- 
opments of work done by others. Typically, Orbits were done 
by Newkirk the year [ was born, Polar plots were done by 
Bishop, and probably others, on machinery resonances in the 
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carly 60's. In fact, Polar plots were quite popular for control 
system work, the area in which I specialized, long before that. 
The Keyphasor™ was used by Stone and Underwood in their 
oil whip work in the mid 40's. The Digital Vector Filter does 
mathematically the fundamental job of convolution that The- 
arle introduced tn his excellent balancing work from the 30's. I 
could continue this list, but perhaps the point 1s clear. 


In the future, inany more changes and improvements in 
machinery and in the observation of this machinery will 
occur. Totally sealed machines are now being built and 
installed. Machines with better efficiency, longer lite, lower 
cost, and higher performance are certain to evolve. The sci- 
ence of the measurement of vibration must also evolve. 


It is worthwhile to reexamine some of the more interesting 
cases In the history of rotor dynamics. This article examines a 
number of early developments and concepts in rotor 
dynamics. The discussion will include the Half-Circle 
Assumption of shatt deflection versus load, bearing lubrica- 
tian pressure versus rotor stability, rotor mass and stability, 
the use of root locus techniques in stability analysis, and 
radial stiffness of fluid-film bearings at zero eccentricity, 
Some of these concepts arc, in the light of our present day 
knowledge, known to be incomplete, Others can be seen to 
have been very forward-looking for their time and can still 
provide useful guidance on potential directions for future 
research. Still others are actively being researched at this 
lime, 


The Half-Circle Assumption of shaft displacement 
versus load 


Qne of the most amusing and unfortunate assumptions 
made in the early days of rotor dynamics is that the plot of 
shaft position versus applied load takes the form of a semicir- 
cle. This semicircle is very crudely (topologically) correct; 


however, It 1s the use of the circle to relate the fluid-film direct 
stiffness to quadrature stiffness that 1s the problem. 
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Figure | shows the load versus displacement assumption. 





Figure | 
Shaft displacement locus inside a cylindrical bearing using the 
Half-Circle Assumption, 


The equation of the assumed half circle is: 





ag 7 ; 

= 4 cy-y= (Vote assumption | 
2 (i) 
where c is the radial clearance, and x and y are the displace- 
ments in the directions of the coordinate axes. Equation (1) 
reduces to: 


xe tye =Cy 
The force equations for this simple system are, for y into x 
(clockwise) rotation (as shown in Figure 1): 


(Note assumption) (2) 


FY, =K + KX (3) 
0=K xK y (4) 
so that, 
Fy 
Y 
K= ——— (5) 
rite + ye 
Fx 
“a x2 + ye i 


Taking the ratio of Equation (6) to Equation (5), 

K x 
a= (7) 
Ka y 








Relating Equation (5) and Equation (6) back to the 
assumed semicircle relationship in Equation (1), these amus- 
ing and unreal relationships are obtained: 





F; 
Ky = (Note assumption) (8) 
+ 
META 
K Se eee (Nate assumption) (9) 
cy 
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The attitude angle, ¥, is actually: 
ix 
+= arctan | — (fd 
4 ¥y 


Therefore, by Equation (7), 


+ = arctan | —— (ii) 





Obviously, the actual relationship between x and y is abso- 
lute as a function of clearance, c, However, with the Half- 
Circle Assumption, while the rigid relationship between 
direct stiffness, Ky, and quadrature stiffness, K,, are now 
forced into the form: 


i S | cy-y? | e 
azy W£ y 





(Note assumption) (12) 








(Note assumption) (£3) 


there is nothing to establish the absolute value of either. It is 
well-known that, for constant rotative speed and fluid viscos- 
ity, the lubricant wedge support stiffness (ic. the quadrature 
stiffness) is essentially constant at low through medium 
eceentricities, so that a plot of Ky and K, may be made, as in 
Figure 2. 

As a typical case, let 


c = .O1 inch (f4) 
Ka = 10,000 Ibin (£3) 
then, using Equation (13), 
y \we 
Ky = 10000 {——— | Ibin (Note assumption) (16) 
0- 
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Bearing Direct Stiffness versus eccentricity ratio comparing typi- 
cal actual behavior to the Half-Circle Assumption. 
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and, using Equation (1), 


x=[(.0l-y}y]'7 in 
Finally, using Equation (8), 
is = Oy kK lb 


(Note assumption) (17) 


(Note assumption) (18) 


The Half-Circle Assumption does not assume that K, 1s 
constant; this additional assumption was made here ta make 
the Ky term look something like its usual nonlinear shape. 

In fact, however, the quadrature stiffness. K, is shown by 
both theary and experiment to be remarkably constant for 
low and medium eccentricity ranges, considering the natural 
nonlinearity and nonsymmetry of typical bearings al constant 
rotative speed, 2. 

Whereas the Half-Circle Assumption can yield a reason- 
able set of Ky and K, terms, the vital error of their relation- 
ship is that the direct dynamic stiffness, Ka. is reasonably 
independent of rotative speed, &. A second error 1s to assume 
that the quadrature dynamic stifiness, K, is directly propor- 
tional to the rotative speed. 

Two secondary errors of the Half-Circle Assumption are 
that (1) the assumed semicircle requires zero direct stiffness 
at the centerline, whereas an externally pressurized bearing 
has high direct stiffness at zero eccentricity, and (2) the fluidic 
inertial effect term may have the exact opposite effect, yield- 
ing a negative direct stiffness at low to medium eccentricities. 
These effects are shown typically m Figure 3. 

The real typical steady displacement versus steady load for 
constant viscosity and various rotative speeds is shown in 
Figure 4. As can be observed, the semicircle is approximately 
correct for just one set of conditions, in the same manner that 
a stopped watch 1s correct twice a day. 


ty 










Dashed Line = Half-Circle Assumption 

EXTERNALLY 

PRESSURIZED BEARING 
Typical 


HIGH FLUIDIC 
INERTIAL EFFECT 
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Figure 3 
Typical shaft locus of static displacement versus static load show- 
ing the effects of external pressure and high fluidic inertial effect. 
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One can conclude that the Half-Circle Assumption can be 
used only for an extremely crude representation of bearing 
behavior. The use of the equation for the semicircular rela- 
tionship in serious rotor dynamics, however, leads to very 
poor results and continued misunderstanding of bearing and 
seal behavior, 


The Lomakin Effect 

Professor A. A. Lomakin wrote a very good paper on Pump 
Stability in 1958. It was his belief, so far as | can interpret, that 
pumps have a quadrature dynamic stiffness which increases 
with the square of speed. Our more recent work shows an 
increase of quadrature dynamic stiffness which ts linear 
(power of one) with speed. Thus, from Lomakin (1958) (our 


Interpretation), 
n 2 
— 1) (Prom Lomakin) (19) 
2 


where Q is the shaft rotative speed, and w is the angular 
velocity of rotor precession. However, we have found (Bently 
& Muszynska (1988)) that: 


Ka = Dja-A n) + D w (20) 


where D is the damping of the fluid bearing, w 1s the angular 
velocity of rotor precession, A is the fluid average circum- 
ferential velocity ratio, and D, is the rotor structural and 
external damping. 

Recent papers on pumps with quotations of that work seem 
to infer that the direct stiffness of a pump increases with the 
square of speed, Surely, there could be such a direct stiffness 


+y 


Low 0 
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Figure 4 
Typical shaft static displacement versus static load showing the 
effect of various rotative speeds, 


_ Orbit 7 


term, especially if the pump pressure increased with the 
square of speed (which it tends to do), and that pressure 
increased the flaw of a seal, bearing, or other rotor-stater 
interface. However, this is, in my belief, not much related to 
Lomakin. 

Itis my guess that it was Lomakin himself who confused the 
issuc. Near the end of the paper, he notes (as translated), 


The frequency of the vibration under apera- 
tonal conditions and the combined action of 
elastic and hydrodynamic forces is: 





P = VPF AP 


[t is my observation that these frequencies do not combine in 
this fashion. Furthermore, it may be seen that his elastic term, 
P, is clearly a direct term; therefore, the hydrodynamic term, 
AP, (as evidenced by the Pythagorean vector addition), 1s a 
tangential (quadrature) term. 


(Note: Lomakin Statement! 


Studies undertaken at Bently Rotor Dynamics Research 
Corporation have shown the existence of a term which we call 
the fluidic inertia effect, because it varies with the square of 
the differential speed, as does Lomakin’s expression. How- 
ever, this effect is negative, not positive, and appears only for 
low and medium eccentricities and has no noticeable effect 
on Stability, either mathematically or in tests. 


Bearing lubrication pressure 

As any technology develops, it occasionally gets off on a 
wrong track, Usually, these incorrect paths are not “bloopers, 
because at the time that they are done they represent good 
work. The problem usually occurs when all followers, in the- 
ory and in practice, continue 1 employ the incorrect line, 
even after it has been clearly replaced by a much better set of 
rules. 

Very typical of this is the use of low oil pressure to feed 
sleeve bearings. It is, even now, repeatedly stated in specifica- 
tions and recommended practices that supply lubricant pres- 
sure be limited, typically to 25 Ibs/in+. The use of very low oil 
pressure helps to assure that the bearing does nor have lubri- 
cant a full 360 degrees around. This tends to make it much 
more difficult for the rotor system to go unstable. This tech- 
nique has been well-known since Harrison (1919) stated that, 
“a fully lubricated bearing is always unstable.” This is, in fact, 
offen correct if nothing is done to correct the situation. 

However, with proper bearing design, a high-pressure 
lubricated bearing is lar superior to a low-pressure fed bear- 
ing and may have much higher efficiency, much better heat 
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transfer, and much better load bearing characteristics. By 
using proper antiswirl and pressurized bearing design, much 
higher (and better controlled) stability characteristics can be 
obtained, 

We have repeatedly shown that the first two cures for fluid- 
Induced rotor system instability sources at bearings are: 

l. Antiswirl 

2. Increased lubrication pressure 

Therefore, this restriction on use of high pressure lubri- 
cants represents a major impediment to progress on rotating 
machinery, 


Rotor mass and stability 

About twenty years ago, it was popular to suggest that, if a 
horizontal rotor system had fluid-induced instability, the 
rotor was too light for its bearing clearance. It has been often 
published that, to cure the instability, more weight must be 
added to the rotor, This rule seems to have quietly (and 
fortunately) disappeared, because adding rotor weight on a 
horizontal machine may have any of three results: 

l. The fluid instability can become worse, 

2. The effect on stability is indifferent. 

3, The rotor system can be stabilized. 

This is because added weight affects both the numerator 
and denominator of the basic natural frequency equation, 








vw, = .|— (21) 


where w, is the natural frequency of the rotor system, Kpr is 
the effective (modal) stiffness of the rotor system, and M „ep 1s 
the effective (modal) mass of the rotor system, 

The effect of the denominator is obvious: if the rotor 
weight, W, goes up, M.e goes up proportionally, and the 
resonant frequency goes down, which is bad for stability. 

The effect of the numerator is also fairly obvious. The 
direct spring coefficient of a bearing is highly nonlinear with 
eccentricity, Increasing steeply at high eccentricity, and the 
direct spring is one of the components of the numerator. 
Thus, when K y (numerator) changes are directly propor- 
tional to Map (denominator) changes, the effect is 
indifferent. 

If Kop grows more rapidly than Mp then stability is 
enhanced, 

However, when adding rotor mass to a horizontal machine, 
very often Mae grows faster than K,,,, so that the fluid 
instability is worsened by the increased rotor weight. 

[tis obvious that the original proponents of this method of 
stability control recognized the distinction between adding 
rotor mass and adding a steady radial sideload. Perhaps later 
authors lost track of this simple distinction. Steady sideload is 
often a good stability enhancement procedure. 


September 1993 


lt is generally beheved that the steady radial sideload pro- 
vided by eravity in horimontal machines is always desirable. 
However, this sideload can lead to several possible problems. 
First, a steady sideload produces a cyclic rolative speed stress 
on the shatt which may contribute to shaft cracking. Second, 
when the shaft has asymmetry and a resonance al or near 
twice rotative speed, each time the stiffer axis of the shafi 
passes the top dead center location of the sideload, the resul- 
lanl snapping action will strongly contribute to shaft fatigue 
and eventual failure, 

Third, much energy may be wasted as heat at the high 
eecentricity spot in the bearing, Lastly, as noted by Kimball 
(1925), Internal Friction Forward Circular Whip Instability is 
encouraged by any detection; therefore, the radial sideload 
cure for instability may also be a cause of instability, Some 
hydrogen and oxygen pumps probably suffer from this 
problem. 


Lund and root locus 

On occasion, proper theory and presentation has been 
developed and only a few pay attention. One example of this is 
the root locus graph done by Jorgen Lund (1974). Lund 
learned from Mel Prohl, and those two are considered by 
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Figure 5 
Root locus plot as presented by Lund (1974). The plot displays 
damped natural frequencies (vertical axis) versus damping expo- 
nent (horizontal axis), The numbers on the curves represent the 
beanng damping coefficients that were vared to analyze the stabil- 
ity of the rotor system he was studying. 
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many, Including myself, to be the major leaders of modern 
rotor dynamics, Figure 5 shows one of the root locus plots 
presented in Lunds paper, 


This is a root locus plot as pioneered by Walt Evans (1954), 
who taught me root locus for control system design. Untor- 
tunately, there are only a few other researchers (Hahn, E. 
Kramer, and T. bwatsubo are among them) who have imple- 
mented this method of stability analysis and followed Lund’s 
lead, 

The root locus is the locus of a family of points, each of 
which represents the roots (cigenvalues) of the characteristic 
equation formed from the set of differential equations of the 
rotor system. These loci can be plotted in a plane representing 
the “real” and “imaginary” parts of the roots. (There is noth- 
ing “imaginary” about these values! They represent real, 
physical machine behavior. Better names for these terms are 
Decay/Growth Rate and Precession Frequency, respectively.) 

The stability threshold (or threshold of instability) is 
reached when the Decay/Growth Rate is equal to zero. Once 
this boundary is crossed, the roots, when based on linear 
differential equations, become meaningless; the machine 
whirl or whip amplitude becomes bounded by nonlinearities 
in the fluid bearing or seal. 

Since rotor system models are based on a large number of 
design paramcters, such as masses, damping coefficients, 
stiffnesses, rolor rotative speed, fluid circumferential velocity 
ratios, and more, any parameter of interest can be varied ina 
systematic way to explore the effect of that parameter on 
stability. 
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Figure 6 
Root locus plot for a typical multi-DOF rotor system model show- 
ing the two lowest modes (both positive and negative roots for 
each mode). The loci represent rotative speed variation from 0 
(marked by the A) to 2000 rac’sec, Only one root poses a threat to 
machine stability because it crosses the vertical axis. 
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Figure 7 
An enlarged view of the important root shown in Figure 6. 
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Figure 8 
Root locus plot showing the effect of antiswirl on the important 
root in Figure 7. Both curves reflect changing rotative speed © (0 to 
2000 rad‘sec) for A = 0.5 and A = 0.1. The instability threshold, 
oth, is about 1000 rad’sec without antiswirl (A = 0.5). With anti- 
swirl (A = 0.1) the instability threshold is increased to greater than 
2000 rad/'sec. 


A typical example of such a root locus plot is shown in 
Figure 6, This graph was created using a multi-DOF, 2-mass 
rotor model which resulted in a sixth-order characteristic 
equation with complex coefficients, and the equation was 
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solved numerically. The rotor rotative speed, 2, was varied 
from Ü to 2000 rad/sec. Only the roots corresponding to the 
lowest two modes are shown. Note that they appear in conju- 
gate pairs at zero rotative speed (at any other rotative speed, 
the roots are no longer conjugate). 

Note also that only one of the four roots shown poses a 
threat to stability and that the nearly conjugate root (reverse 
precession) root moves toward greater stability. This is almost 
always the case; the first mode forward precession root is 
usually the only one that poses a stability problem. Figure 7 
shows an enlarged view of the important root that was shown 
im Figure 6, 

Root locus techniques can be used to study the effects of 
changing design variables on the stability and damping factor 
of a machine at the design stage and later to evaluate fixes to 
stability problems discovered during machine operation. 


Figure & shows the results of an antiswirl study on the root 
shown in Figure 7. The rotor speed was again varied from 0 to 
2000 rad/sec. The fluid average circumferential velocity ratio, 
A, was varied from 0.5 to 0.1. The dramatic improvement in 
stability is evident in the figure. The machine, which was going 
unstable at a rotative speed of around 1000 rad/sec, is now 
slable at a rotative speed of over 2000 rad/sec, 

The root locus is a powerful tool for stability control and 
will be used heavily in the future. Bently Rotor Dynamics 
Research Corporation 1s actively researching the application 
Of rool locus techniques to the solution of rotor dynamics 
stability problems, 


Poritsky and radial stiffness 


Another case of a man ahead of his time is that of H. H, 
Poritsky. He did the basic equations of stability very well for 
his time. He, as well as everyone else, used 1/2 as the lubricat- 
ing fluid average circumferential velocity ratio which, 
although tog restrictive, was not fundamentally wrong. How- 
ever, his excellent work (Poritsky 1953) was heavily criticized 
because he included a radial spring at zero eccentricity. 
Apparently, other researchers claimed he was incorrect on 
this point, thereby disregarding his entire contribution, It is 
true that it is possible to have zero radial spring stiffness at 
zero eccentricity. I have built such bearings many times, and, 
of course, because the resulting mechanical resonance is at 
zero speed, they allow the machine to go into oil whirl at zero 
rotative speed. 

However, most ordinary bearings do have at least a small 
radial stiffness at zero eccentricity, or alternatively, the steady 
state rotor position 1s slightly eccentric, or the shaft or seals 
supply a direct spring component. Thus, Poritsky’s work was 
actually extremely useful, 


Conclusion 


History can be useful and interesting, but I do not believe in 
steering a boat by watching the wake. This is a partial review 
of things that went wrong, things that went right, and some 
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ideas about what will occur in the 
Future. This study of “land mines” was 
made as an adjunct of another study 
which will be the genesis and anteced- 
ents of the Quadrature Dynamic Stift- 
ness term published by Dr. Agnes 
Muszyoska and myself, 


Kq = (w-AnyD. (22) 


I am pleased that it has a genesis, 
because tt has been extremely difficult 
to get this paradigm shill inta general 
acceptance, 

Lack of knowledge and lack of mea- 
surement of these dynamic stiffness 
parameters represent a major handicap 
to the rotating machinery industry, This 
is because, without them, it is imposs- 
ible to measure and specify the perfor- 
mance characteristics necessary to 
prevent instability. 

This article shows some of the history 
ol rotor dynamics in a form, hopefully, 
that will be useful to others in their pur- 
suit of future improvements to our 
knowledge and control of machinery 
behavior. W 

Fditors Note: The December T993 
Orbit will include an article written by 
Donald Bently on root focis. 
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The ratio of 2X to 1X vibration - 
A shaft crack detection myth 


by Malcolm Werner 
Corporate Manager, MES 
Bently) Nevada Corporation 


or many years, a myth has cir- 
culated in industry that 
relates the presence of and/or 


growth of a shaft crack to the ratio of 


two times synchronous vibration 
amplitude (2X) to a synchronous 
vibration amplitude {1X}: 


. 2X Amplitude 
Crack Ratio = = e 
a 1X Amplitude 





The myth states that a crack is pre- 
sent or a crack is growing when this 
ratio is increasing. If the ratio is con- 
stant or if there is no 2X, then there 
isnt a shaft crack, according to this 
myth, 

Bently Rotor Dynamics Research 
Corporation (BRDRC) has been 
“keeping score” on shaft cracks for 


more than ten years, Based on this 
“score card” of actual shaft crack 
cases, If you follow the crack 
ratio formula, you will be exactly 
wrong more than 75% of the time. 


BRDRC, Bently Nevada's 
Machinery Diagnostics Services, and 
our customers, have been using a 
field-proven shaft crack diagnostic 
methodology for many years. We 
have found that, in more than 75% of 
documented cases, the kev indicator 
of shaft cracks is change in 1X ampli- 
tude or 1X phase lag angle. While 2X 
vibration is helpful in crack diagnos- 
tics, In more than 254 of cases there 
was no 2X vibration. 


In addition, there are other 
sources for 2X vibration, such as the 
effect of a steady sideload into a non- 
linear stiffness bearing or seal when 
there is some 1X bow. Therefore, the 
crack ratio can not only fail to predict 
a crack, it can give false indication of 
a crack. 

The “score card” summarizes the 
key shaft crack indicators, based on 
experience. W 


Vibration score card for shaft cracks 


Changing parameter 


1X Amplitude 


1X Phase 


Prevailing indicator i in 1 shaft crack cases 


75% of cases. Nearly always increasing 
| amplitude 


| Often seen with 1X amplitude, probably rare by 
| itself, Can be leading, lagging or constant phase, 


ss, Either 24 Amplitude (usually 


increasing) or Amplitude and Phase changes, 


2X Phase 


2X Amplitude 25% of c 





depending on rotor system and crack properties. 
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